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Abstract 

A class of inhomogeneous nonlinear Schrodinger equations (NLS), claiming to be novel integrable 
systems with rich properties continues appearing in PhysRev and PRL. All such equations are shown 
to be not new but equivalent to the standard NLS, which trivially explains their integrability features. 

PACS no: 02.30.1k , 04.20.Jb , 05.45.Yv , 02.30.Jr 

Time and again various forms of inhomogeneous nonlinear Schrodinger equations (IHNLS) along 
with their discrete variants are appearing as central result mostly in the pages of Phys. Rev, and 
PRL CLU76,RB87,PRA91,Kon93,PRL00,PRL05,PRL07, which are either suspected to be integrable 
due to the finding of particular analytic or stable computer solutions, or assumed to be only Painleve 
integrable arxiv08, or else claimed to be completely new integrable systems. Apparently the solution of 
such integrable systems needs generalization of the inverse scattering method (ISM), in which the usual 
isospectral approach involving only constant spectral parameter A has to be extended to nonisospectral 
flow with time-dependent A(i). Moreover certain features of the soliton solutions of such inhomogeneous 
NLS, like the changing of the solitonic amplitude, shape and velocity with time were thought to be 
new and surprising discovery. 

We show here that all these IHNLS , though completely integrable are not new or independent 
integrable systems, and in fact are equivalent to the standard homogeneous NLS, linked through simple 
gauge, scaling and coordinate transformations. The standard NLS is a well known integrable system 
with known Lax pair, soliton solutions and usual isospectral ISM nls,ALM. As we see below, a simple 
time-dependent gauge transformation of the standard isospectral system with constant A can create the 
illusion of having complicated nonisospectrality. Similarly, a time-dependent scaling of the standard 
NLS field Q — > q = p(t)Q would naturally lead the constant soliton amplitude to a time-dependent 
one. In the same way a trivial coordinate transformation x — > X = p(t)x would change the usual 
constant velocity v of the NLS soliton to a time-variable quantity v(t) = and the invariant shape 

of the standard soliton with constant extension T = ^ to a time-dependent one with variable extension 
r(i) = (see Fig la a,b). Therefore all the rich integrability properties of the IHNLS, observed 
in earlier papers, including more exotic and seemingly surprising features like nonisospectral flow, 
appearance of shape changing and accelerating soliton etc. can be trivially explained from the time- 
dependent transformations of these IHNLS from the standard NLS and the corresponding explicit result 
, namely the Lax pair, N-soliton solutions, infinite conserved quantities etc. for the inhomogeneous NLS 
models can be derived easily from their well known counterparts in the homogeneous NLS case through 
the same transformations nls. 
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Let's start from a recent version of IHNLS PRL07, which is generic in some sense: 

iQt + \DQ XX + R\Q\ 2 Q - {2ax + fl 2 x 2 )Q = 0, (1) 

where the nonautonomous coefficients D(t),R(t) of the dispersive and the nonlinear terms are arbitrary 
functions of t and the other time-dependent functions are 

a(t) = s t p, n 2 (t) = \{9 t - l -D6\ where p = ^, 9 = ^, (2) 

s(t) being another arbitrary function. It is easy to see that a time-dependent scaling of the field can 
change the coefficient of the nonlinear term in ([1]) and at the same time generate an additional term 
from iQt, while a change in phase of the field involving x 2 would yield extra terms from Q xx . As a 
result transforming Q —* q = yfpe l i x Q. we can rewrite IHNLS (pQ) into another form 

Ht + yfc + D\q\ 2 q - (2ax + iD9)q - iD9xq x = 0. (3) 

In PRL07 Eq ([T]) was declared to be a new discovery and as a proof of its integrability a Lax pair 
associated with Eq ([3|) was presented, which we rewrite here in a compact and convenient form by 
introducing a matrix U^°\q) = yfa{qo + — q*a~) as 

U(X(t)) = -iX(t)a 3 + U^(q), 

V(X(t)) = DV (X(t)) -iaxa 3 + D9xU(X(t))), (4) 

where 

V (X(t)) = -i\ 2 (t)a 3 + X(t)U {0) + \^{U^ - (C/(°)) 2 ). (5) 

We can check from the above Lax pair that the flatness condition Ut — V x + [U, V] =0 yields the IHNLS 
([3]) under the constraint X(t)t = a + D9X(t). Using relations ([2]) one can resolve this constraint to get 
X(t) = p(t)(X + s(t)), which was given in PRL07. We now establish the equivalence between the Lax 
pair U(X(t)), V(X(t)) dH[5|) for the IHNLS and the well known Lax pair U n i s (X), V^ S (A) of the standard 
NLS n]s, showing explicitly that the nonisospectral X(t) is convertible to constant spectral parameter 
A through simple transformations. For this it is interesting to notice first, that the structure of the 
NLS Lax pair is hidden already in the expression of the IHNLS Lax pair as U(X(t) = A) = U n i s (X) 
and Vo(X(t) = A) = V^ S (A). Therefore the aim should be to remove the t-dependence from X(t) = 
p(t)(X + s(t)) by absorbing the arbitrary functions p(t) and s(t) in step by step manner. Note that 
the Lax pair U, V, as evident from the associated linear problem J|<I> = = V$, correspond 

to infinitesimal generators in the x and the t direction, respectively and therefore a simple coordinate 
change (x,t) -> (x = p(t)x, t = t) resulting Jl = ±j|, ^ = JL + ptxjj* = ^ + D °x-§^, would yield 
U(x,t) = pU(x,t) and V(x,t) = V(x,t) + D9xU(x,t). Therefore using such a transformation and 
comparing with @, we can easily remove the p(t) factor from X(t) in U(x,t), which however would 
scale the field as q — > 2 and at the same time eliminate from the transformed V(x,t) the nonstandard 
term D9xll(x,t) appearing in V( x, t) For the removal of additive term p(t)s(t) from X(t), present 

~ ' 3 

in U(x,t), one can perform a gauge transformation $ — > $ = <7<J> with g = e tpsa , taking the Lax pair 
to a gauge equivalent pair 

U = g x g- x + gUg-\ V = g t g- x + gVg' 1 . (6) 

One notices that though the above transformations are enough to remove explicit t dependence from 
U due to its linear dependence on X(t), the removal of t from V(X(t)) becomes a bit involved due to 
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the nonlinear entry of \ 2 (t) and X(t)U^ in it, which bring in more time-dependent terms like 2Dp 2 s 
and Dp 2 s 2 . These extra terms however can be exactly compensated for by extending slightly the above 
coordinate and gauge transformations by introducing additional functions f(t), f(t) and choosing them 
as ft = 2Dp 2 s and ft = Dp 2 s 2 . The multiplicative factor Dp 2 appearing in all terms in V(X(t)) can be 
absorbed easily by a further coordinate change t — ► T = Dp 2 t. Therefore taking the above arguments 
into account one finally solves the problem completely through the following three steps of simple 
transformations: 

i)Coordinate transformation : (x,t) — > (X,T), X = p(t)x + f(t), T = Dp 2 t, with ft = 2Dp 2 s, (7) 
it) Gauge transformation §6§, where g = ^(pW^M /(*))«* w j t h f t = Dp 2 s 2 , (8) 
iii)Field transformation : q — > ip, where ip = ~TjX*l e 2 '^^*^* 4 "^*^ . (9) 
The above transformations would take ([!]) directly to the standard NLS Lax pair 

C/ nis (A) = -iAd 3 + U (0 \ where U (0) = ipa + - ^V", 

V nls (X) = ^a. + XU^+'-asiuP-iU^) 2 ). (10) 

which proves the equivalence of the Lax pair Q for the IHNLS ([3]) and the Lax pair (|1U|) associated 
with the standard NLS: 

#T + -^XX + = 0, (11) 

obtained as the flatness condition of (j 1 U j) . One can also check that under the change of independent and 
dependent variables ([7]) and ([9]) the inhomogeneous NLS ([3]) is transformed directly to the homogeneous 
NLS (|TTjh 

Therefore we remark that the inhomogeneous NLS ([I]) and ([3]) are equivalent to the homogeneous 
NLS (jlip . a well known integrable system. The corresponding Lax pairs ((4|) and (jlOl) are also gauge 
equivalent to each other, which therefore trivially explains the complete integrability of the inho- 
mogeneous NLS. All signatures of the complete integrability like the Lax pair, N-soliton solutions, 
infinite conserved quantities etc. for these HNLS can be obtained easily from the corresponding well 
known expressions for the NLS system (1111 [TUl) by inverting the set of transformations (I7|8|9j) as 
(X,T) — > (x, t), g — > g" 1 , ifi —> q. As a result, explicit t-dependence obviously enters in the Lax 
operators as well as in the amplitude, phase and the x-dependence of the field q{x) of the IHNLS 
system, resulting the spectral parameter A — > X(t) and making the constant amplitude A, extension 
r and velocity V of the soliton to become t-dependent. Fig. 1 demonstrates this situation, show- 
ing that the NLS soliton (module) a) = ^sech^, £ = ±(X - VT) goes to IHNLS soliton b): 
\q\ = A(t)seche, i = ^(x - v(t)), where A(t) = Ap(t), v(t) = DVp(t)t - j§,T(t) = ^ under 
the transformations inverse to (|7l8l9p . Therefore even though IHNLS soliton (Fig. lb) looks rather 
exotic and quite different from the standard NLS soliton (Fig. la), these solutions are related simply 
by coordinate and scale transformations and belong to equivalent integrable systems. 

It is worth mentioning that, though in all earlier papers only 1-soliton of the IHNLS was consid- 
ered, one can easily derive the exact N-soliton for the IHNLS, thanks to its complete integrability, by 
exploiting again its equivalence with the integrable NLS, i.e. by simply mapping the known N-soliton 
of the standard NLS through the same transformations ([7119]). 

By redefining the field further: Q — > b(t)Qe ia ^ with arbitrary functions a(t),b(t), we can generate 
more inhomogeneous terms in ([TJ resulting a more general form of IHNLS 

iQt + jQxx + R\Q\ 2 Q-(2ax + n 2 x 2 + a + i~f)Q = 0, j(t) = j, (12) 
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equivalent naturally to the integrable NLS. The IHNLS (|12h was found to be the maximum inhomo- 
geneous NLS system which can pass the Painleve integrability criteria PLA87. A recently proposed 
IHNLS arxiv08, which is simply a particular case of (|12[) at a = and a = 0, is therefore also equivalent 
to the standard NLS and hence, contrary to the assumption in arxiv08 that the system is only Painleve 
integrable and not completely integrable, the equivalence with NLS assures the complete integrability, 
including the existence of infinite conserved quantities, N-soliton solutions etc. for this IHNLS arxiv08 

We now look into other forms of integrable HNLS appeared earlier in Phys. Rev. RB87,PRA91,Kon93 
and PRL CLU76,PRL00,PRL05 and show their equivalence to the standard NLS, similar to as found 
above. The simplest form of inhomogeneity to the NLS: 2xq was proposed in CLy76, which is clearly 
a particular case of (J3J) with a = 1,0 = 0, ensured by the choice R = D = 1, s = t, proving thus its 
equivalence with the NLS. 

A more general IHNLS with F(x)Q was considered in RB87 and shown finally that integrability 
restricts the choice only upto F(x) = a + ax + fix 2 , which is consistent with the general integrable 
IHNLS (|12p , shown to be equivalent to the standard NLS (jlip . However for constructing such integrable 
IHNLS, as shown here, x-dependent spectral parameter considered in RB87 is not needed and similarly 
the restriction on function h(t) appearing in X(t), found by the author apparently as a condition for 
the integrability, actually does not appear allowing the function to be arbitrary, as shown here. 

In PRLOO a variant of IHNLS was considered, which was suspected to be integrable through com- 
puter simulation. It is easy to see however, that this IHNLS can be obtained as a particular case 
from (fT2|) at a = 0, s = 1, a = 0, 7 = 0, £1 = 0, but with nontrivial R(t),D(t),^(t) obeying certain 
constraints. Similarly IHNLS proposed in PRL05 can be seen to be derivable from (|12j) as a particular 
case with D = l,R(t) = goe ct ,s = I, a = 0,7 = 0, giving a = 0, but £1 = —c 2 . Therefore both 
these inhomogeneous NLS PRL00iPRL05 are equivalent to the standard NLS and hence completely 
integrable. 

Some integrable discrete versions of IHNLS, namely inhomogeneous Ablowitz-Ladik models (ALM) 
were proposed in PRA91,Kon93, containing in addition to the standard ALM ALM an inhomogeneous 
term nojip n , with u = 1 PRA91 or u(t) as an arbitrary function Kon93 . We find that in spite of 
the discrete case a similar reasoning found here holds true and the proposed inhomogeneous ALM can 
be shown to be gauge equivalent to the standard ALM ALM, under discrete gauge transformation: 
U n = gn+iUng' 1 , V n = gnVng' 1 + gnQn 1 ! witn 9n = e~ mF ^ a and redefinition of the field as q n -> 
iftn = Qn^ 2n+1 ^T(t), where Tt(t) = uj(t) is an arbitrary function as found in Kon93. 

Based on the above result we therefore conclude that the general inhomogeneous NLS, if integrable, 
should be of the form (|12p . Other forms of integrable IHNLS are only its particular cases. However 
all these inhomogeneous NLS are neither new nor independent integrable systems, but are equivalent 
to the standard homogeneous NLS, from which all their integrable structures like Lax pair, N-soliton 
solutions, infinite number of commuting conserved quantities etc. can be obtained easily through simple 
mapping. The time-dependent soliton amplitude, shape and velocity as well as the nonisospectral flow 
in these inhomogeneous NLS are just an artifact of the time-dependent coordinate, gauge and field 
transformations, needed to get these systems from the standard NLS. Therefore before proposing any 
new integrable inhomogeneous NLS the authors should check whether it can be linked in any way to 
the general integrable IHNLS (I12p . whose equivalence with the well known NLS we have proved here. 
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Figure 1: Exact soliton solutions (module) for the integrable a) homogeneous NLS (jlip and b) inho- 
mogeneous NLS © with A = 2.0, V = 0.5. and the particular choice R = .008t 2 ,D = l,s = 0. In 
spite of the significant differences between the appearance and dynamics of these two solutions they 
are related by simple transformations ([Till]) and belong to equivalent integrable systems. 
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